Reionization - A probe for the stellar population and the physics of the early universe. 
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We calculate the reionization history for different models of the stellar population and explore the 
effects of primordial magnetic fields, dark matter decay and dark matter annihilation on reionization. 
We find that stellar populations based on a Scalo-type initial mass function for Population II stars 
can be ruled out as sole sources for reionization, unless star formation efficiencies of more than 
10% or very high photon escape fractions from the parental halo are adopted. When considering 
primordial magnetic fields, we find that the additional heat injection from ambipolar diffusion and 
decaying MHD turbulence has significant impact on the thermal evolution and the ionization history 
of the post-recombination universe and on structure formation. The magnetic Jeans mass changes 
the typical mass scale of the star forming halos, and depending on the adopted stellar model we 
derive upper limits to the magnetic field strength between 0.7 and 5nG (comoving). For dark 
matter annihilation, we find an upper limit to the thermally averaged mass-weighted cross section 
of {av)/moM < 10"'^'' cm'^/s/eV. For dark matter decay, our calculations yield a lower limit to the 
lifetime of dark matter particles of r > 3 x 10^"^ s. These limits are in agreement with constraints 
from recombination and provide an independent confirmation at a much later epoch. 
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I. INTRODUCTION 



Over the last decades, our understanding of cos- 
mological reionization has improved considerably. 
Observations of high-redshift quasars clearly indicate 
that reionization must end around ~ 6 [l|, and the 
WMAP 5-year-measurement finds a reionization optical 
depth of T = 0.087 ± 0.017 0, [3], yielding clear evidence 
that reionization started significantly earlier and is 
thus a continuous process. Reionization through mini- 
quasars would produce a soft-X-ray background that 
is significantly higher than the observed background, 
and can thus be ruled out [4|, Q . Simulations of stellar 
reionization further show that the process is highly 
inhomogeneous and based on the growth and merging of 
ionized bubbles around the first stellar sources 6, 7, 8]. 
Understanding the essential physical processes, it is 
possible to build semi-analytic models that describe 
stellar reionization, which can be tested over a wide 
parameter space [I, El, [H [H [H, [11. Additional 
constraints on the cosmic star formation rate are now 
available from gamma-ray burst observations [l^|. The 
reionization framework thus provides an increasingly 
reliable test for the stellar population during reionization 
and can be used to constrain global physical conditions 
in the early universe. 



|18l | and [l9| indicate that gas especially in more massive 
systems can fragment because of dips in the equation 
of state, which may lead to the formation of a stellar 
cluster. It was shown that cooling in previously ionized 
gas is enhanced and leads to typical stellar masses of 
~ 10 Mq no, [m. It was further suggested that the 
presence of weak magnetic fields is sufficient to lead to a 
more present-day like mode of^star formation, resulting 
in considerably lower masses |23i | . Recently, a new phase 
of stellar evolution was suggested in which the stars 
would be powered by dark matter annihilation instead 
of nuclear fusion |24l |. Various follow-up works have 
explored such a scenarios in more detail. Studies by 
locco [231, Freese et al. Q explored the main-sequence 
and pre-main-sequence phase of dark stars, and other 
works calculated the effect of dark matter capture by 
off-scattering from baryons [13, HI HI, [ll . While many 
predictions are still model-dependent, it has often been 
suggested that such stars may have typical masses of 
800 Mq, giving rise to a very bright main-sequence 
phase, or may have much longer lifetimes due to 
modifications in the stellar evolution. In the end, all 
these suggestions must face the constraint that the 
stellar population must be able to provide the correct 
reionization optical depth. 



Regarding the stellar population, the works [la] and 
[l7| suggest that the first stars were top-heavy with a 
peak in the IMF at around 100 Mq. On the contrary. 
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Complications may arise through the presence of ad- 
ditional physics that are often not considered in stan- 
dard reionization calculations and simulations on the first 
stars. Such possibilities include the presence of primor- 
dial magnetic fields, dark matter decay and dark matter 
annihilation. Magnetic fields have been observed on all 
scales in the universe, and recently, it has been demon- 
strated that they were present already in high-redshift 



galaxies (3l|. They are found in the interstell ar g as as 
well as in the intergalactic medium [32, [Sj, [3J], but 
their origin is still unclear. There is a viable possibil- 
ity that these fields have a primordial origin 35, 36]. 
So far, the most stringent constraints on the strength of 
these putative fields come from the measurements of the 
cosmic microwave background radiation (CMBR) and 
from big-bang nucleosynthesis (BBN) calculations. A 
homogeneous magnetic field would produce temperature 
anisotropics in the CMBR [33| whose maximal ampli- 
tudes are limited by the COBE satellite measurements 
which in turn limit the field strength, as measured to- 
day, to Bo < 3.5 X 10~^ G [13]. The presence of any pri- 
mordial field would also alter the CMBR power spec- 
trum by changing the characteristic v eloc ities. With 
the sensitivity of the PLANCK satellite [13lj one should 
be able to detect fields with present day strength of 
So > 5 X 10~^ G i3g]. So far, measurements by the 
WMAP satellite [132] are compatible with the absence 
of primordial magnetic fields. 

Strong magnetic fields in the early universe can also 
change the abundance of relic ''He and other light ele- 
ments during the big bang nucleosynthesis [33, |40] • To 
comply with observational limits on light element abun- 
dances these primordial fields must not exceed 10^^ G at 
the time when the universe was T = 5 x 10^ K which 
corresponds to a present day field i3o^3x 10~^ G[4Q]. 

Effects of primordial magnetic fields have already been 
considered in |41| , finding that density perturbations can 
be enhanced by the Lorentz force from tangled magnetic 
fields. Both the evolution of perturbations in the pres- 
ence of magnetic fields as well as their effect on the ther- 
modynamics via ambipolar diffusion heating and decay- 
ing MHD turbulence was considered in [12] . Recent cal- 
culations of [43] show that the enhancement of structure 
due to magnetic fields is pronounced at about 5x 10^ Mq, 
but becomes less effective on larger mass scales and ap- 
pears as a subdominant contribution on the scale of the 
magnetic Jeans mass. Consequences for 21 cm observa- 
tions have been explored as well [ij, Sa] • In this work, 
we examine the consequences for reionization in more 
detail and calculate the backreaction on structure forma- 
tion according to the work of [6|. The WMAP 5 year 
data [S] have measured the Thomson scattering optical 
depth from reionization and allow to constrain different 
reionization scenarios in the early universe. 

The nature of dark matter is still unclear, and conse- 
quences of various particle physics like massive neutrinos 
or axion decay have been explored early [e. g. l46l. l47l. [48| . 
Observational progress allowed to refine these studies and 
to explore such scenarios in the framework of ACDM cos- 
mology m [HO, [U, m, [H, H [El. Recently, the possi- 
bility was discussed to detect such effects using future 
21 cm telescopes (53, and it was shown that the fraction 
of the energy absorbed into the IGM can be calculated 
in detail for specific models of dark matter decay and 
annihilation [57|. Indeed, secondary ionization through 
the decay / annihilation products can provide a way to 



ionize the IGM which is independent of the stellar contri- 
bution |58j], whereas the additional heat input increases 
the Jeans mass and delays the formation of the first struc- 
tures. 

In this work, we show how reionization constrains the 
properties of the stellar population as well as some addi- 
tional heat sources like primordial magnetic fields, dark 
matter annihilation and decay. In ijlll we present our 
model for stellar reionization, which considers the IGM 
as a two-phase-medium consisting of ionized bubbles and 
overall neutral gas. Based on the thermal evolution, we 
self-consistently determine the minimal mass scale of ha- 
los which can collapse. In mill we explain our treatment 
of primordial magnetic fields and show how they mod- 
ify the thermal evolution. The treatment of dark matter 
annihilation and decay, as well as some implications for 
dark stars, are discussed in mVl The optical depth for 
different models is given in ^V\ Further discussion and 
outlook is given in WIl 



II. REIONIZATION IN THE EARLY UNIVERSE 

The thermal and ionization history of the IGM be- 
tween recombination and the end of reionization is deter- 
mined by a number of different processes. At redshifts 
z > 300, Compton scattering of CMB photons couples 
the gas temperature T to the CMB temperature Trad- 
At lower redshifts, this coupling is less efhcient and the 
gas decouples from the CMB due to adiabatic expansion. 
This standard scenario can be altered if additional energy 
injection mechanisms are present. Additional heat input 
can lead to an earlier redshift of decoupling, whereas an 
increase in the ionized fraction, for instance due to sec- 
ondary ionization from the decay or annihilation prod- 
ucts of dark matter, tends to make Compton cooling 
more efficient. Once structure formation sets in, the 
thermal and ionization history is further influenced from 
X-rays produced in star forming regions and UV pho- 
tons that escape from the first galaxies. We modified 
the RECFAST code [H and included all these feedback 
processes as described below, modeling the IGM as a 
two-phase medium of ionized and partially-ionized gas. 
In this picture, the fully-ionized gas refers to the gas in 
the HII regions of the first luminous sources, while the 
partially ionized gas describes gas that is not yet affected 
by UV feedback. 



A. The RECFAST code 

In the absence of additional energy injection mecha- 
nisms, recombination in the early universe and freeze- 
out of electrons was calculated with unprecedent accu- 
racy, following the detailed level populations of hundreds 
of energy levels for H, He and He'*' and self-consistently 
calculating the radiation field[60|. They developed the 
RECFAST code |l33l ] , a simplified version of the multi- 



level calculations, based on an effective three-level model 
for the hydrogen atom. RECFAST is capable of fully 
reproducing the results of the more detailed calculation 
[5g |. Both the detailed calculation and the RECFAST 
code were recently updated and improved in ^61] . For the 
partially-ionized gas, we modify the equations describing 
the thermal and ionization history in the following way: 
The equation for the temperature evolution is given as 
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where Lhcat is the new heating term (see ^III Al mil Bl 
and mi Cp , Lcooi the new cooling term including Lyman a 
cooling, bremsstrahlung cooling and recombination cool- 
ing using cooling functions of [62], ax is the Thom- 
son scattering cross section, ur the Stefan-Boltzmann 
radiation constant, rUe the electron mass, c the speed 
of light, ks Boltzmann's constant, n the total number 
density, Xe ~ Uc/nn the electron fraction per hydrogen 
atom, H{z) is the Hubble factor and /hc is the num- 
ber ratio of He and H nuclei, which can be obtained as 
/hb = ^/4(1 — Yp) from the mass fraction Yp of He with 
respect to the total baryonic mass. The evolution of the 
ionized fraction of hydrogen, Xp, is given as 
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In this equation, nu is the number density of hydro- 
gen atoms and ions, hp Planck's constant, fcjon is the 
collisional- ionization coefficient [63| , /ion describes ion- 
ization from X-rays (see ijIIC[) . and the parametrized 
case B recombination coefficient for atomic hydrogen an 
is given by 



au^F X 10 
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l + ct'^ 



cm^ s ^ 



(3) 



with a = 4.309, b = -0.6166, c = 0.6703, d = 0.5300 and 
t = r/10^ K, which is a fit given in [6J] to the coefficient 
of [6a|- This coefficient takes into account that direct 
recombination into the ground state does not lead to a 
net increase of neutral hydrogen atoms, since the photon 
emitted in the recombination process can ionize other 
hydrogen atoms in the neighbourhood. The fudge factor 
F = 1.14 serves to speed up recombination and is de- 
termined from comparison with the multilevel-code. We 



further introduce the clumping factor C{z) 



l)/{r 



to take into account the increase in the recombination 
rate in structures of increased density at low redshifts. 
We use the fit formula 



C{z) = 27.466exp(-0.114z + 0.001328^2 



(4) 



obtained from simulations of '661 ^^ redshifts z < 40 and 
set C{z) = 1 at higher redshifts. 

The photoionization coefffcient /3h is calculated from 
detailed balance at high redshifts as described in 59, 60]. 
At lower redshifts, however, radiative recombination is 
no longer balanced by photoionization in the presence 
of additional energy injection mechanisms like ambipo- 
lar diffusion heating. Once the ionized fraction drops 
below 98%, we thus calculate the photoionization co- 
efficient directly from the photoionization cross section 
given in [63] ■ The frequency i'H,2p corresponds to the 
Lyman-a transition from the 2p state to the Is state of 
the hydrogen atom. The two-photon transition between 
the states 2s and Is is close to Lyman-a. Consequently 
we use the same frequency for both processes. Finally, 
Ah — 8.22458 s"-'^ is the two-photon rate for the tran- 
sition 2s-ls according to [68| and K^ = A^ 2p/[^^F[{z)] 
the cosmological redshifting of Lyman a photons. The 
additional terms for Eq. ^ and ^ that are needed to 
describe the effects of dark matter annihilation and decay 
will be discussed in mV\ 

The fully-ionized component in the HII regions is de- 
scribed with the volume filling factor Qh+ that denotes 
the volume fraction of ionized hydrogen bubbles. It is 
needed to compute the reionization optical depth and 
takes UV feedback into account. It evolves as 



dQii+ _ QB_+C{z)nc^}i+aA duph/dz 



dz 



H{z){l + z) 
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(5) 



d 0, [O, [H, [11 [11. The UV photon production 
rate duph/dz will be described in miC\ in more detail, 
no,H+ denotes the mean electron number density in fully- 
ionized regions, nn the neutral hydrogen density in re- 
gions that are still unaffected from UV photons and 
a A = 4.2 X IQ-^^Tmax/lO'^ K]-°-^ cmVs is the case A 
recombination coefficient [63 which we have chosen here, 
as recombinations will preferably occur in dense regions 
where the recombination photons are unlikely to escape 
into the IGM [70j. It is evaluated at the temperature 
with T„iax — max(10''' K, T) to account for the fact that 
the ionized regions should be heated at least to lO'* K. If 
heating via ambipolar diffusion and decaying MHD tur- 
bulence increased the gas temperature above this thresh- 
old, it is evaluated at the gas temperature instead. To 
compare the models with the observational constraints 
from WMAP, we calculate the Thomson scattering opti- 
cal depth of free electrons, given as 
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where nH,o is the comoving number density of ionized 
and neutral hydrogen, f2A and flm the usual cosmological 
density parameters, Hq the Hubble constant and Zrc is 
the redshift where reionization starts, which we define 
as the point where the effective ionized fraction X(,ff = 
Qh+ + (1 — Qii+)xp becomes larger than 6.5%. This 
point has been chosen after a comparison of the actual 



TE cross calibration spectra with those assuming sudden 
reionization, but the results are not very sensitive to this 
choice. In summary, we adopt the following picture: 

• The partially-ionized gas unaffected by UV feed- 
back is modeled with Eqs ([T]) , ([3]) that describe the 
evolution of its temperature and ionization degree. 

• The fully-ionized gas is described by the volume- 
fiUing factor Q^+ calculated from Eq. ([5]). Its tem- 
perature is given as the maximum of 10"* K and the 
temperature of the partially-ionized gas. 

• Recombination both in the partially- and fully- 
ionized gas is enhanced by the clumping factor 
given in Eq. ^. 

• For the partially- ionized component, recombina- 
tion is described using the case B recombination 
coefficient, appropriate for low-density gas far from 
the ionizing sources. 

• For the fully-ionized gas, we use the case A recom- 
bination coefficient, assuming that most recombi- 
nations take place in high-density regions near the 
ionizing sources. 

The implementation of different feedback mechanisms is 
described below in more detail. 



B. The generalized filtering mass 

The universe becomes reionized due to stellar feedback. 
We assume here that the star formation rate (SFR) is 
proportional to the change in the fraction of collapsed 
halos /coll. As shown in [7i|, /coii is given as 



/coll = erfc 
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and where rrimin is the minimum mass of haloes that 
are able to collapse at a given redshift. Here, ii{x) — 
(sinx — xcosx)/x'^ is the spherical Bessel function, R is 
related to the halo mass Mh via Mh — 47rpi?'^/3, p is 
the mean density, 5c ~ 1.69/L'(z) the linearized density 
threshold for collapse in the spherical top-hat model and 
D{z) the linear growth factor. In the absence of magnetic 
field, m,„i,i is determined by the filtering mass p, |72| . 
As discussed in the introduction, tangled magnetic fields 
can potentially create more small scale structure via the 
Lorentz force. The calculations by 43| show that this 
effect is pronounced in the minihalo regime. While they 
adopted a constant minimal collapse mass of l(fh^^ M© 
which is independent of the magnetic field, we use the 
framework of [6l| to take into account the change in the 



mass scale of halos that can form stars. Indeed, we 
find that the mass scale is changed significantly (see Fig. 
1(d)), such that the additional small-scale structure from 



tangled magnetic fields does not contribute to star for- 
mation. We introduce the magnetic Jeans mass, [42] 
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The filtering mass in the presence of magnetic fields is 
then given as 



M 



2/3 
F,B 



da'Ml'^a') 



a 



1/2- 



(11) 



where a = 1/(1 + z) 



is the scale factor and Mg — 



m.ax{Mj,Mj). The minimum halo mass to consider is 
given as mmin = max(Mi;'_B,7TT,cooi), where TOcqoI denotes 
the minimal halo mass for which the baryons can effi- 
ciently cool after collapse. We adopt here the fiducial 
value mcpcL = 10^ Mq 17311 . Similar estimates have been 
given in [7J], [ll| and [73. For the subsequent analysis, 
we furthermore assume that only a certain fraction f^ of 
the collapsing halo mass turns into stars. 



C. Stellar feedback 

As X-ray photons have long mean free paths, they can 
play an important role in the ionization and heating of 
the gas. Assuming that the local correlation between the 
SFR and the X-ray luminosity (from 0.2 — 10 keV) holds 
up to a renormalization factor fx [see \WL VTA U&L VM. ISOJl . 
the X-ray heating function Lx is given as 



Lx = 3.4 X IO^Va- 



SFR 



1 Mq yr- 



erg s 



(12) 



X-ray emission has two major sources, inverse-Compton 
scattering of CMB photons with relativistic electrons ac- 
celerated in supernovae, and high-mass X-ray binaries. 
The former may play an increasingly important role at 
high redshifts, as the CMB photons are more energetic 
at high redshifts [84I. In the early universe the factor 
fx ^ 0.5 if 10^^ erg are released per 100 Mq, correspond- 
ing to an overall efficiency of 5% |81| . The abundance of 
X-ray binaries depends on metallicity and the stellar ini- 
tial mass function, and could be especially large if very 
massive Pop. Ill stars dominate 80]. In the model pre- 
sented here, we adopt fx ~ 0.5 as a generic value and 
assume that the uncertainty can be ascribed to the star 
formation efficiency f^. In general. X-rays lose their en- 
ergy through three channels. A fraction fx,h goes into 
heating, a fraction /x,ion mto ionization, and a fraction 



/x,coii into excitation. These parameters should not be 
confused with /ion introduced in Eq. ^ and /cou in- 
troduced in Eg ([7|). We calculate them using the fit 
formulae in [118J . We can thus write 
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The contribution from the X-rays to the term /ion in Eq. 
Q is then given as 



/ic 
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(14) 



UV photons from stellar sources are in general absorbed 
locally, only a fraction fcsc that manages to escape from 
the first galaxies can contribute to the ionization of the 
IGM. Again, the production of ionizing photons is asso- 
ciated with the star formation rate |76|, yielding a con- 
tribution 



where 



duph/dz ^ df coll 
Uh dz 



Z, — ^Hc/* /csc-'»ionj 
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(16) 



with Ayic = 4/(4 - 31^) = 1.22 and iVion is the number 
of ionizing photons per stellar baryon. 



D. Models for the stellar population 

The nature of the first stellar sources is still under dis- 
cussion, and even the question whether the primordial 
IMF is top-heavy or closer to the locally measured IMF 
is not solved. While for instance tlfil and ^^ found a 
top-heavy IMF using adaptive-mesh refinement (AMR) 
or smoothed-particle hydrodynamics (SPH) simulations, 
it was shown in [l8| that primordial and low-mctallicity 
gas can fragment if the evolution of the gas is followed 
further after the formation of the first clump, due to a 
dip in the equation of state. Based on similar arguments, 
|ig | argued that even a small metallicity fraction can lead 
to fragmentation in the first protogalaxies. It was further 
suggested that magnetic fields can have a crucial influ- 
ence on the primordial IMF [23| . Given these uncertain- 
ties, we discuss different models for the stellar popula- 
tion. Dark stars have also been suggested as some of the 
first luminous sources. We will discuss this possibility in 
more detail in ijIVCl 

In models A and B, we assume that reionization is 
solely due to metal-free massive Pop III stars. This sit- 
uation can only be examined from a theoretical point of 
view. Most investigations of UV feedback from high-mass 
zero-metallicity stars indicate that ionizing photons can 
easily escape the star- forming halo and drive large HI I 
regions into the IGM. This suggests that the escape frac- 
tion in high-redshift galaxies could be very high, of order 



A 


III 


40,000 1 


B 


III 


40,000 1 (Tvi,. < lO"* K) 
0.1 (Tvir > 10'' K) 


C 


III/II 


10,000 1 (Tvir < lO"* K) 
0.1 (Tvir > lO-'K) 


D 


II 


4,000 0.06 



TABLE I: Summary of adopted stellar models. The first col- 
umn gives the model name. The second column indicates the 
stellar populations that contribute to reionization. The third 
column gives the number of ionizing photons per baryon used 
in Eq p6|l . The fourth column lists the adopted escape frac- 
tions. Model A is a highly extreme case in which we assume 
all ionizing photons come from massive Pop III stars and can 
escape the star-forming halo. In model B photons can escape 
efficiently only from halos less massive than a corresponding 
virial temperature of 10* K, while for higher-mass halos, this 
fraction is reduced to 10%. In model C we assume that Pop 
III as well as Pop II stars contribute to reionization, and we 
adopt the same escape probabilities as in B. Model D is an- 
other extreme case which assumes that all ionizing photons 
come from low-mass Pop II stars, with escape fraction 6%. 



unity, being much higher than in the present-day universe 
[for instance [SJ, Hj, [S^l . Detailed numerical simulations 
of ,8G] show that indeed the shock bounding HII regions 
of massive Pop. Ill stars can easily photo-evaporate the 
minihalo and lead to an escape fraction of one. However, 
the situation is not fully clear. There are other studies 
that suggest that the escape fraction could remain small 
[Sq . depending on the mass of the collapsing halo. This 
is why we adopt two different models. Model A is an 
extreme case, in which we assume that the escape frac- 
tion is 100% independent of the halo mass. The more 
realistic case is probably model B, where we distinguish 
whether the virial temperature corresponding to the gen- 
eralized filtering mass is smaller or larger than 10"* K. If 
it is smaller, we still assume that the halo is easily pho- 
toevaporated. If it is larger, we assume that most of the 
mass is kept within the halo and adopt an escape frac- 
tion /esc = 0.1. The virial mass corresponding to 10^ K 

is given as Af^ = 5 x IO^Mq {J^\ [tI, [s^ ■ For model 

A and B, we adopt a total number of ioniz ing photons 
per baryon of A^ion = 40, 000, following [12^ . 

The heavy elements produced by the very first stars 
will gradually enrich the IGM. It is very likely that some 
contribution to reionization comes from low-metallicity 
Pop II stars as well. As these stars are expected to have 
lower masses than Pop HI stars, we introduce model C, 
which corresponds to a stellar population of intermedi- 
ate mass. It has A'ion = 10, 000 and an escape fraction 
according to model B, suggesting that such stars may 
photo-evaporate only low-mass halos. 

Finally, we consider in model D the extreme case of 
very rapid chemical enrichment and assume cosmic reion- 



ization is entirely driven by low-mass Pop II stars. We 
assume a stellar mass distribution similar to the local 
IMF in the Galactic halo [H, S, ^ with a metal- 
licity of 1/20 of the solar value. This corresponds to 
-^ion ^ 4, 000. We furthermore adopt escape fractions 
that are typical for local star forming galaxies. Many 
upper limits and a few detections have been found ob- 
servationally [93, 94, 95, 96, 97], suggesting /osc "^ 0.06. 
A detection from Lyman-break galaxies at z ^ 3 implied 
a higher fraction of 10% [92], while more recent observa- 
tions place upper li mits of 5 — 10% or cl aim detections 
at even lower levels [lOi, [lol [lol [lOi]. As the rele- 
vant rates scale only linear in /esc , it seems reasonable to 
adopt /esc = 0.06 for our Pop II model. 

As it is unlikely that cosmic reionization occurs in- 
stantaneously with the onset of Pop III star formation 
nor that chemical enrichment is so rapid that all ionizing 
photons come from metal-enriched Pop II stars, we adopt 
the intermediate scenario C with a mixed population as 
our fiducial model. All global cosmological parameters 
are chosen according to the WMAP 5 year data Q. 



III. THE EFFECT OF MAGNETIC FIELDS ON 
THE IGM 

As discussed in [42, [93], magnetic fields can signifi- 
cantly alter the thermal evolution of the IGM via am- 
bipolar diffusion heating and decaying MHD turbulence, 
which can have a significant influence on the filtering 
mass scale. We explain our treatment of these heating 
terms in this section. 



A. Ambipolar diffusion heating 



and He , which is a good approximation in the absence 
of shocks. They are given by 



(o-w)h+,Hc = (1.424 + 7.438 X 10~^T 
- 6.734 X lO^^T^) X 10 



(19) 
(20) 



-9 cm3 s"i. 



As the power spectrum for the magnetic field is unknown 
and eq. (J17p cannot be solved exactly, we adopt a simple 
and intuitive approach to estimate the integral for a given 
average magnetic field B with coherence length L. The 
heating rate can then be evaluated as 
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IQn'^jplpt L"^ 



(21) 



The coherence length L is in principle a free parame- 
ter that depends on the generation mechanism of the 
magnetic field. It is constrained through the fact that 
tangled magnetic fields are strongly damped by radia- 
tive viscosity in the pre-recombination universe on scales 
smaller than the Alfven damping scale k^l^^ given by 
[T07 l[T08l . [l09t 



234 Mpc" 



n^h^^'^^ 
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10-9 G 






1/4 



0.02 



h 
OJ 



1/4 



(22) 



where Bq = B/{1 + z)^ denotes the comoving magnetic 
field. In fact, we expect fiuctuations to be present on all 
scales. As the heating term goes as L~^, we thus adopt 
the minimal value L — fcmax/(l + 2)- 



The presence of magnetic fields introduces two differ- 
ent contributions to the heating rate, one coming from 
ambipolar diffusion and one resulting from the decay of 
MHD turbulence. 

In the first case, the contribution can be calculated as 
112,191: 



L 



anibi 



16n^Jpf,p^ 



V X B] X B 



(17) 



Here, p„, pi and pi, are the mass densities of neutral 
hydrogen, ionized hydrogen and all baryons. The ion- 
neutral coupling coefficient for primordial gas is given as 

poirioii 
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\nH{(Jv) 



H+,H 



|nHe(0-u)H+,He 



ruH Jn-H + 4nHe] 



(18) 



where nue is the number density of He and mu the mass 
of the hydrogen atom. Collisions with electrons are ne- 
glected here, as their contribution is suppressed by a fac- 
tor m^/mH- We adopt the z ero drift velocity momentum 
transfer coefficients of [1061 for collisions of H+ with H 



B. Decaying MHD turbulence 

For decaying MHD turbulence, we adopt the prescrip- 
tion of [42], 



B^itfia [\n{l + U/U)fH{t) 
8tt 2 [\n{l + td/U)+\n{t/U)]'^^ 



(23) 



where t is the cosmological time at rcdshift z, td is the 
dynamical timescale, U the time where decay starts, i. e. 
after the recombination epoch when velocity perturba- 
tions are no longer damped by the large radiative viscos- 
ity, Zi is the corresponding redshift. For a power spec- 
trum of the magnetic field strength with power-law in- 
dex aj_the parameter d is given as a = 2{a + 3)/(a -I- 5) 
(llOl IlllL III2I Ill3| ]. In the generic case, we expect the 
power spectrum of the magnetic field to have a maximum 
at the scale of the coherence length, and the heat input 
by MHD decay should be determin ed from the positiv e 
slope corresponding to larger scales [ll2l . Ill3l . Ill4 illSl ] . 
We thus adopt a = 3 for the calculation. We estimate the 
dynamical timescale as td = L/va, where va = B / ^A-Kph 



is the Alven velocity and pb the baryon mass density. 
The evolution of the magnetic field as a function of red- 
shift can be determined from the magnetic field energy 
Eb = -B^/Stt, which evolves as H 



dEE 
dt 



— —4:H{t)EB — iambi " L, 



decay- 



(24) 



C. The evolution of the IGM 

While the redshift of reionization of course depends 
on the model for the stellar population, the mechanism 
which delays reionization in the presence of magnetic 
fields is always the same. We illustrate this for our fid- 
ucal model C adopting a star formation efficiency /* of 
1%, but point out that the general discussion is valid also 
for the other stellar population models. We use the cos- 
mological parameters of WMAP 5 [2^ . Fig. 1 (a) shows 



the thermal evolution of the IGM for different magnetic 
field strengths. In the absence of magnetic fields, the gas 
temperature follows the temperature of radiation until 
z ^ 200, where Compton scattering becomes inefficient. 
Afterwards, the gas cools adiabatically until it is reheated 
during reionization. In the presence of magnetic fields, 
gas and radiation decouple earlier due to the ambipo- 
lar diffusion heating and decaying MHD turbulence and 
stays at higher temperatures, which prevents collapse in 
low-mass halos and thus delays reionization. The addi- 
tional heat increases the ionized fraction at early times 
(see Fig. 1(b) I, while the redshift where the IGM be- 



comes fully ionized is delayed, due to the increased gen- 
eralized filtering mass (see Fig. 1(c) I. Collisional ion- 



ization introduces a natural temperature plateau of the 
order ^ 10000 K, as any further temperature increase 
will have a backreaction on the ionized fraction of the 
gas, and thus make ambipolar diffusion less efficient. 



IV. IMPLICATIONS FROM THE DARK 
SECTOR 

A large variety of different particle physics models has 
been sug geste d to explain the dark matter content of the 
universe J116l | . Many of these models predict some inter- 
actions in the dark sector and include some form of dark 
matter decay and annihilation. The consequences of such 
scenarios for the thermal evolution have been discussed 
in detail for instance in [5y, [53|. In this work, we discuss 
the implications of such scenarios on reionization. As 
we have seen for the case of primordial magnetic fields, 
the additional heat input increases the filtering mass and 
thus increases the minimal halo mass in which the first 
luminous objects can form, making stellar reionization 
less effective. On the other hand, once heating through 
dark matter decay or annihilation is effective, secondary 
ionization will likely be effective as well and thus increase 
the Thomson scattering optical depth. As a consequence 
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FIG. 1: Evolution of different quantities as a function of red- 
shift, a) The effective gas temperature, b) The effective ion- 
ized fraction, c) The filtering mass, d) The comoving mag- 
netic field strength. 
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of dark matter annihilation models, a new phase of stel- 
lar evolution was suggested in [2J, [2^ , in which stars are 
not powered by nuclear fusion, but by annihilating dark 
matter within them. These stars could reach masses of 
up to 10'^ Mq [2y|. Once dark matter annihilation be- 
comes ineffective and such stars enter the main sequence 
phase, they would thus be extremely bright sources of UV 
photons. The consequences of such a phase are discussed 
below. Further models taking into account a dark matter 
capturing phase have been suggested as well [23, [23, [30] ■ 
As they are highly parameter-dependent, we study them 
in a separate work [ll7[. 



for which we use the fitting formulae of 118^. As the de- 
cay rate scales with the dark matter density squared, it 
is most efficient at early times and may thus modify the 
recombination history, which allows to place upper limits 
on the dark matter annihilation. As shown in [119] , the 
WMAP 1-year-data yield an upper limit of 



(av) < 2.2 X 10 



-29 en? s 



-1 „-i r mDM \ 

Jabs \ 7i,r„T/ i 



V MeV I ' 



(28) 



where jahs corresponds to the energy fraction actually ab- 
sorbed into the IGM and is of order 1 for particle masses 
in the MeV range [57| . 



Dark Matter annihilation 



B. Dark Matter decay 



The energy of the annihilation products can be de- 
posited into heating, ionization and coUisional excitation. 
The latter is quickly radiated away, but may contribute 
to the build-up of a Lyman a background that helps to 
couple the spin temperature of hydrogen to the gas tem- 
perature, which may be relevant for 21 cm observations 
[5a |. For given models of dark matter, it is possible to 
work out the detailed absorbed fractions as a function of 
redshift [57|. In this work, however, we adopt the more 
generic approach of [76|. For definiteness, they assume 
that the dark matter particles annihilate to high-energy 
photons, which allows to determine the energy fractions 
going into heat and ionization by the fitting formulae of 
[ll8J. As they show, heating by dark matter annihila- 
tion is accompanied also by a significant increase in the 
electron fraction due to secondary ionization. For other 
models of dark matter annihilation, the results can be 
rescaled by appropriate absorption efficiencies. To Eq. 
P]) describing the temperature evolution, we thus add a 
term 



8T 



m'm-pC 



3mkB{l + z)H{z) 



ixxh 



(25) 



according to [56|, where 771 = 1 + /hc + Xi, rj2 = 1 + 4/ho, 
/ho is the helium fraction by number, Xh is the fraction 
of energy going into heating [see 118] and S,x the effective 
baryon-normalized "lifetime" , given as 



e 



X 



rriDM 



{av){l + zf 



n 



DM 



(26) 



Here, ^dm and fif, denote the cosmological parameters 
for dark and baryonic matter, p^ is the critical density at 
redshift zero, muM is the dark matter particle mass and 
{av) the velocity-averaged cross section. In the same way, 
we add a term to Eq. ([3]) that describes the evolution of 
the ionized fraction: 



5xi fmpC^ 



(27) 



where -Bion — 13.6 eV is the hydrogen ionization thresh- 
old and Xi th6 fraction of the energy going into ionization. 



In a similar way, dark matter decay can alter the ther- 
mal and ionization history of the universe. The lifetimes 
of dark matter particles considered here will be consider- 
ably larger than the age of the universe to ensure that the 
abundance of dark matter does not change significantly. 
Still, the conversion of particle mass to thermal energy 
can have an important impact on the thermal evolution of 
the universe [50, [53]. In this Hmit, the decay rate is con- 
stant over time and the effect at lower redshifts is more 
pronounced compared to the case of dark matter anni- 
hilation. To calculate the heating and ionization rate, 
the same formalism can be employed as for dark matter 
annihilation, but we adopt the baryon-normalized decay 
rate 
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DM 
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^btx ' 



(29) 



where tx is the lifetime of the dark matter particle. 
Based on the modified recombination and reionization 
histories in the presence of dark matter decay, it was 
shown that Il20ll 



JxJO,bs 

tx 



< 2.4 X 10^ 



-25 



(30) 



where f^ is the fraction of the particle mass-energy re- 
leased through decay and fabs the fraction of the released 
energy deposited into the IGM. If data of large-scale 
surveys are included, the constraint can be slightly im- 
proved, yielding 



tx 



< 1.7 X 10^2^ 



Dark stars 



(31) 



Dark stars have been suggested as a possible conse- 
quence of dark matter annihilation in the early universe 
|24l \23l . [2q . If indeed the first stars form in the cen- 
tral peaks of the dark matter distribution in the early 
minihaloes, they might be powered by dark matter anni- 
hilation instead of nuclear fusion. Such stars could have 



up to 10'^ Mq and surface temperatures in the range 
of 3000 — 10000 K, i. e. significantly colder than con- 
ventional Pop. Ill stars. However, this dark phase will 
gradually come to an end as the dark matter in the cusp 
annihilates, so the star may contract and finally reach 
a main sequence phase. A linear stability analysis and 
ID simulations including hydrodyna mics and radiation 
indicates the stability of such objects [l2ll |. In this case, 
one can expect a similar number of photons per stellar 
baryon as for conventio nal P op. Ill stars, perhaps even 
higher by a factor of 2 [l22j |. As it is crucial for these 
objects to form on the peak of the dark matter density, 
we expect them to form only in low-mass halos with a 
virial temperature below 10^ K, as the so-called atomic 
cooling halos are very turbulent f73'| and are more likely 
to form a stellar cluster instead of a single massive star 
|18| . In addition, such massive halos are more likely to 
have accreted material from previous metal enrichment 
|123i] . which may cause fragmentation as well 19]. In pre- 
viously ionized regions, the thermodynamics of collapse 
are significantly altered and simulations generally find 
lower-mass stars [20l |2l| , which might be an important 
limitation for dark stars as well. 
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FIG. 2: The reionization optical depth for Pop. Ill stars 
(model A) in the presence of primordial magnetic fields. The 
contour fines are equally spaced around Tro ~ 0.87 with 
At = 0.017, corresponding to the difi'erent a-errors of the 
measurement. Magnetic fields higher than 20 nG are clearly 
excluded by the optical depth alone. The transition between 
stellar reionization and coUisional ionization occurs at about 
5 nG, providing a more stringent limit if metal enrichment is 
required. 



V. CONSTRAINING THE PARAMETER SPACE 
WITH WMAP 5 

In this section, we show the results for the different 
stellar populations and in the context of different addi- 
tional physics. 



A. Stellar reionization with and without 
primordial magnetic fields 

We first concentrate on the effect of primordial mag- 
netic fields and compare the results of different models 
for the stellar population. We run our reionization model 
for a range of different star formation efficiencies and 
primordial magnetic fields and obtain the corresponding 
optical depth. To constrain the parameter space, some 
assumptions need to be made. The first is obvious: The 
calculated optical depth must agree with the optical 
depth measured by WMAP 5 at least within 3a. In 
addition, we require that the star formation efficiency 
may not be unreasonably high. For Pop. Ill stars 
with ~ 100 Mq forming in halos of -- 10^ M© (dark 
matter), star formation efhciencies of the order 0.1% 
are expected. It seems thus reasonable to reject star 
formation efficiencies higher by one order of magnitude, 
i. e. of the order 1%. Further requirements are that 
reionization must be complete by redshift 6, and that 
there should be a stellar contribution to reionization, 
so that the universe becomes metal-enriched. Now we 
discuss the constraints resulting from these criteria. 
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FIG. 3: The reionization optical depth for Pop. Ill stars 
(model B) in the presence of primordial magnetic fields. Mag- 
netic fields higher than 20 nG are clearly excluded by the op- 
tical depth alone. The transition between stellar reionization 
and coUisional ionization occurs at about 5 nG. We can fur- 
ther exclude magnetic fields between 2 and 5 nG within 3-a, 
as too high star formation efficiencies would be required to 
obtained the measured optical depth. 



For all stellar models, we can exclude magnetic fields 
larger than 20 nG, based on the optical depth alone 
(see Fig. [3]to[S|). a transition occurs at about 5 nG: 
For lower magnetic fields, reionization is still mostly 
due to stellar radiation, whereas for higher magnetic 
fields, the IGM is heated to such high temperatures that 
coUisional ionization becomes very efficient and increases 
the optical depth. In this regime, structure formation is 
strongly suppressed, as the magnetic Jeans mass scales 
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FIG. 4: The reionization optical depth for a mixed popula- 
tion (model C) in the presence of primordial magnetic fields. 
Magnetic fields higher than 20 nG are clearly excluded by the 
optical depth alone. The transition between stellar reioniza- 
tion and coUisional ionization occurs at about 5 nG. In addi- 
tion, magnetic fields between 0.7 and 5 nG can be excluded, 
as they would require unreasonably high star formation effi- 
ciencies for the measured optical depth. 
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FIG. 5: The reionization optical depth for Pop. II stars 
(model D) in the presence of primordial magnetic fields. For 
magnetic fields lower than 5 nG, stellar reionization domi- 
nates, but unreasonably high star formation efficiencies would 
be required to reconcile the measured optical depth. For 
stronger magnetic fields, coUisional reionization dominates, 
which can be excluded if metal-enrichment is imposed. Mag- 
netic fields larger than 20 nG can be excluded from the optical 
depth alone. 



with Bq and dominates over the thermal Jeans mass. 
Independently of the steUar model, such a regime is 
always found between 5 and 20 nG. Therefore, structure 
formation is considerably impeded such that metals do 
not form and the universe is not fully ionized by redshift 
6. We can thus exclude magnetic field strength above 5 
nG independent of the stellar model. 



FIG. 6: The reionization optical depth for Pop. Ill stars 
(model B) in the presence of dark matter annihilation. At 
{av)/mDM ~ 10~^^ cm^/s/eV, secondary ionization of the 
annihilation products starts to dominate over stellar reioniza- 
tion. Values higher than 3 x lO^'^"' cm^/s/eV are ruled out 
by the WMAP 5-year data. 



If we further require that the star formation effi- 
ciency must be lower than 1%, the results become 
model-dependent. In model A, which is an extreme 
case where very high escape fractions are assumed even 
for very massive halos, no further constraint on the 
magnetic field is possible on the 3a level. Still, even in 
this extreme case, one can reject fields larger than 2.5 
nG on the la level. The other extreme case, model D 
with low escape fractions and a low number of ionizing 
photons per stellar baryon, can be rejected completely, 
as it always requires unphysically high star formation 
efficiencies. 



For more realistic cases, essentially model B and C, we 
find that a critical magnetic field strength exists above 
which a very high star formation efficiency is needed 
to get into the 3a interval around the measured optical 
depth. For model B, the critical value is 2 nG, for model 
C, it is at 0.7 nG. As it seems likely that a transition 
to less massive Pop. II stars occurs during reionization, 
model C might indeed be the most realistic case and pro- 
vide an upper limit of 0.7 nG. 



B. Dark reionization scenarios 

For the case of dark matter annihilation and decay, we 
will show only the results for the Pop. Ill star model B. 
Regarding the competition between stellar reionization 
and secondary ionization from high-energetic annihila- 
tion / decay products, this corresponds to a rather con- 
servative choice if one assumes that stellar reionization 
should do the main contribution. Such a stellar popula- 
tion could be interpreted either as conventional Pop. Ill 
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FIG. 7: The reionization optical depth for Pop. Ill stars 
(model B) in the presence of dark matter decay. For lifetimes 
lower than 3 x 10^** s, secondary ionization of the annihi- 
lation products starts to dominate over stellar reionization. 



Lifetimes lower than 3 x 10 
5-year data. 



s are ruled out by the WMAP 



stars or as so-called dark stars in a main sequence phase. 
We find that dark matter annihilation starts to become 
important for {(jv)/mDM ~ 10~^^ cm^/s/eV and com- 
pletely dominates the stellar contribution for higher val- 
ues (see Fig. [5]). Values higher than 3 x lO^'^'^ cm'^/s/eV 
are inconsistent with the WMAP 5-year data. The re- 
sults for dark matter decay are given in Fig. [71 Decay 
starts to become important for lifetimes below 3 x 10^* 
s. However, lifetimes smaller than 3 x 10^^ s are incom- 
patible with WMAP 5. As pointed out recently, con- 
straints on dark matter annihilation can be significantly 
improved taking into account the d ark m atter clumping 



factor C{z) = (Pdm(^))/(/Odm(^))^ [l2j|. However, this 
quantity is highly uncertain and may vary by six orders of 
magnitude [l25|. In the framework of light dark matter, 
we recently showed that is is constrained to be smaller 
than 10^ at redshift zero |l26l |. yielding an enhancement 
of the annihilation rate by a factor of 10 or more at red- 
shift 20. Other dark matter models Il25j | even yield a 
clumping factor that drops down to 1 at these redshifts, 
such that our calculation altogether provides a conserva- 
tive and firm upper limit. 

From Fig. [SI we can further see that very massive dark 
stars in a main sequence phase, which should have higher 
star formation efficienes that are an order of magnitude 
higher than conventional Pop. Ill stars, are ruled out 
by WMAP 5, as the models require rather low star for- 
mation efficiencies to reproduce the Thomson scattering 
optical depth. Of course, this conclusion holds only if one 
assumes that reionization is completely due to dark stars, 
which seems rather unlikely (see discussion in i jIVC|) . and 
we will discuss some alternatives in fjVT] 



VI. CONCLUSIONS AND OUTLOOK 

We have calculated the reionization optical depth for 
different stellar models in the presence of primordial mag- 
netic fields as well as dark matter annihilation and de- 
cay. The results indicate which star formation efficiencies 
are required in the presence of some additional heating 
mechanism for a given stellar model. Considering differ- 
ent stellar models and primordial magnetic fields, we find 
the following results: 

1. Independent of the model for the stellar population, 
we can securely exclude primordial magnetic fields 
larger than 5 nG. 

2. For the most realistic case with a mixed stellar pop- 
ulation (model C), we even find an upper limit of 
0.7 nG, as higher magnetic fields would require star 
formation efficiencies larger than 1%, which is un- 
realistic. Similar results are found for model B, 
assuming reionization completely due to Pop. Ill 
stars. 

3. Reionization only due to population II stars (model 
D) is ruled out completely. 

For dark matter, we found the following results: 

1. Dark matter annihilation provides noticeable con- 
tributions to the reionization optical depth only for 
thermally averaged mass-weighted cross sections 
{(Jv)/mDM > 10-33cm3/s/eV. 

2. Parameters (o'u)/?71dj\/ > 3 x 10^^'^ cm^/s/eV can 
be ruled out by 3a on the basis of WMAP 5 year 
data. 

3. Dark matter decay becomes important for the 
reionization optical depth for lifetimes below 3 x 

1024 S. 

4. Dark matter lifetimes below 3 x lO^'^ s are ruled out 
by 3cr. 

These results are in agreement with conservative con- 
straints obtained from the gamma-ray background [127| . 
We further showed that reionization can not be due to 
~ 1000 M© dark stars alone, as the corresponding optical 
depth would be significantly too high. One might won- 
der whether heating from dark matter annihilation might 
help to significantly delay stellar reionization, in order to 
reconcile this model with observations. However, as can 
be seen in Fig. [SI this is more than compensated by 
the effects of secondary ionization, once that dark mat- 
ter annihilation starts to have a significant influen ce on 
the IGM. In case collider experiments like the LHC jl34| 
or other dark matter detection experiments [135] find ev- 
idence for a self-annihilating dark matter candidate, this 
can be seen as evidence for a rapid transition towards 
a different mode of star formation, or a problem in our 
understanding of dark stars. To reconcile dark stars with 
observations, the following scenarios seem feasible: 
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1. Dark stars with 1000 M© may be extremely rare 
objects, and their actual mass scale is closer to 
the mass scale for typical Pop. Ill stars. In such 
a case, reionization could not distinguish between 
dark stars and conventional Pop. Ill stars. 

2. The transition to lower-mass stars occurs very 
rapidly, such that dark stars only contribute to the 
very early phase of reionization. Reasons for that 
might be chemical, radiative as well as mechanical 
feedback (see also the discussion in i jIVC|) . In fact, 
even a double-reionization scenario might be con- 
ceivable, in which very massive dark stars ionize 
the universe at high redshifts. Due to chemical and 
radiative feedback, formation of such stars might 
be suppressed and the universe might become neu- 
tral again, until the formation of less massive stars 
becomes efficient enough to reionize the universe. 

3. Dark stars do not reach a main-sequence phase, but 
are disrupted earlier by some non-linear instability. 
Such an instability would be constrained by the fact 
that it neither appears in a linear stability analy- 
sis nore in ID simulations. Violent explosions of 
dark stars might even be considered as a source for 
Gamma-Ray bursts. 

While too definite conclusions on the existence of dark 
stars are not yet possible, it is at least indicated that the 
possibilities mentioned above should be explored in more 
detailed, and a better understanding of the properties 
of the dark stars. A better understanding of their evo- 
lution after the dark phase will certainly help to better 
understand their possible role during reionization. 

The constraints derived here are independent of other 
works that essentially rely on the physics of recombina- 
tion to derive upper limits on additional physics [e. g. 
[22l . Ill9l . Il20l |. In particular, magnetic fields can evolve 
dynamically and their field strength may thus change be- 
tween these epochs. It has recently been suggested that 
the Biermann Battery effect creates magnetic fields in the 
presence of an electron pressure gradient [l28i] . We thus 
need to probe magnetic fields at different epochs. As we 
will show in a separate work [45|, upcoming 21 cm mea- 



surements will allow to probe the thermal history before 
and during reionization in great detail, and may allow to 
detect primordial magnetic fields of the order of 0.1 nG. 

Additional ways of probing the reionization history and 
the dark ages exist as well: Scattering of CMB-photons 
in fine-structure lines of heavy elements may lead to a 
frequency-dependent GMB power spectrum and may al- 
low t o me asure metal-abundances as a function of red- 
shift [l29|. Before reionization, molecules may form in 
the IGM and introduce further frequency-dependent fea- 
tures in the GMB |l30l | . Such features are likely enhanced 
in the presence of either primordial magnetic fields or 
dark matter decay / annihilation, as the increased elec- 
tron fraction catalyses the formation of molecules, and 
the additional heat input leads to a departure of the level 
populations from the radiation temperature. 

Further improvements are expected from the upcoming 
measurement of Planck, that will measure the reioniza- 
tion optical depth with unprecedent accuracy and thus 
allow to strengthen the constraints obtained here on the 
stellar populations and additional physics. In addition, 
a more accurate determination of the cosmological pa- 
rameters will remove further uncertainties in the present 
models of reionization. 
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